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Abstract 

The density profile of cool core of intracluster gas is investigated, for a cluster of galaxies that is initially 
in the virial equilibrium state and then undergoes radiative cooling. The initial gas profile is derived under 
the assumption that the gas is hydrostatic within the dark matter potential presented by so-called NFW 
or King model and has a polytropic profile. The contribution from masses of gas and galaxies to the 
potential is ignored compared to the dark matter in the calculation. The temperature and density profiles 
of gas in its quasi-hydrostatic cooling phase, which is expected to last for ~Gyr, is then calculated for 
different initial gas profiles. It is found that in the quasi-hydrostatic cooling phase, while the temperature 
decreases to be about one-third, the density increases by a factor of 4-6 at the cluster center in comparison 
with their initial polytropic values, though the profiles over the core depend on the dark matter potential. 
Hence, the core radius in the quasi hydrostatic cooling gas appears smaller than the initial polytropic one. 
Wc compare the density profile of the cool core with observations to find that while the initial density 
is around the upper bounds of large-core (> 100 kpc) clusters, likely most relaxed but cooling is not yet 
significant, the central density under quasi- hydrostatic cooling falls between the mid- and high- values of 
small-core (< 100 kpc) or cool-core clusters. It is also found for the quasi-hydrostatic cooling gas that the 
entropy profile roughly agrees with the best-fit model for the ACCEPT cluster sample with a low central 
entropy, and the pressure gradient in the inner core is close to that of the REXCESS sample. X-ray surface 
brightness calculated for the quasi-hydrostatic cooling gas is well represented by the conventional double 
/3-model, giving a physical basis of applying the double /3-model to cool core clusters. 

Key words: galaxies: clusters: general - cooling flows - galaxies: intergalactic medium - X-rays: 
galaxies: clusters 



1. Introduction 

Clusters of galaxies offer excellent laboratories for ther- 
mal evolution of cosmic baryons: the objects are filled 
with hot (T ~ 10^ — 10^ K), highly- ionized intracluster 
gas and the baryon fraction in massive clusters is close to 
the cosmic value. Since the timescale of radiative cooling 
at the center of regular clusters is estimated to be shorter 
than the Hubble time, their core regions are thought to 
be affected by cooling. According to statistical studies 
of clusters using X-ray observations, the fraction of clus- 
ters having a compact, cool core (often termed CC clus- 
ters), is roughly 50% (e.g., O'Hara ct al. 2006; Ota et 
al. 2006; Chen et al. 2007; Santos et al. 2008; Cavagnolo 
et al. 2009; Hudson et al. 2010). 

Based on earlier works on peaked surface-brightness 
that cooling core clusters exhibit, it was suggested that 
global cooling-flow would occur unless some heating pro- 
cess balances with radiative cooling (e.g., Fabian 1994). 
The similarity and smoothness of cooling profiles indi- 
cate need for continuous, distributed heat source (for re- 
view, e.g., Peterson & Fabian 2006). On the other hand, 
a possible mild inflow with a quasi-hydrostatic balance 
is proposed to account for observed temperature profiles 
of cooling cores (Masai & Kitayama 2004). Later, using 



a hydrodynamics code, Akahori & Masai (2006) demon- 
strate that such a cooling phase lasts for ^ Gyr and then 
breaks into the global cooling-flow. They discuss the ori- 
gin of observed two-peaked core size distribution (Ota & 
Mitsuda 2002; Ota & Mitsuda 2004) in view of the thermal 
evolution. 

While intracluster gas has often been studied on the 
basis of the conventional isothermal /3-model, a deviation 
therefrom has been commonly observed at the center of 
cool core clusters: they exhibit a systematically higher 
central density while the profiles are fairly universal out- 
side 0.1r,5oo ~ 100 kpc (e.g., Neumann & Arnaud 1999). 
To better reproduce those X-ray observations, some au- 
thors have introduced empirical models such as double-/? 
model (e.g., Jones & Forman 1984) and modified /3 model 
with a cool density cusp at the center and steepening at 
a large radius (Vikhlinin et al. 2006). This model gave 
a good fit to the high-resolution Chandra data including 
the core emission (see also Bulbul et al. 2010). Although 
those models are given in analytic forms and therefore sim- 
ply utilized for observational studies, physical background 
that they are based on is not always clear. 

In the present paper, wc intend to obtain a fiducial den- 
sity profile of cooling gas under quasi-hydrostatic balance. 
Masai & Kitayama (2004) give a temperature profile for 



2 



N. Ota, K. Onzuka, & K. Masai 



the initially isothermal gas, but no density one. We cal- 
culate the temperature of cooling gas in the same manner 
but consider more generally that the gas is initially hy- 
drostatic with an polytropic profile in the NFW (Navarro 
et al. 1997) or King dark matter potential. In section 2, 
our model calculation is described, and pressure and en- 
tropy as well as density of quasi-hydrostatic cooling gas 
are presented. In section 3, we compare those calculated 
quantities with observations, and also examine application 
of widely used /?-model to X-ray surface brightness calcu- 
lated for quasi-hydrostatic cooling gas. Wc then discuss 
implications on the ICM thermal evolution in the cluster 
core regions. Since the profiles of gas density etc. of cool 
cores are of our interest, we present the quantities nor- 
malized to their central values unless specified otherwise. 

2. Calculations 

2.1. Initial hydrostatic gas distribution 

Wc consider a spherically symmetric galaxy cluster ini- 
tially in the virial equilibrium. We ignore the mass of the 
gas and galaxies compared to the dark matter, which is 
thought to occupy typically more than 80% of the cluster 
mass. Thus, the gas distributes so that its pressure bal- 
ances with the local gravitational potential formed by the 
dark matter. For non-cool core (NCC), i.e, for the initial 
state before cooling becomes significant, we assume that 
the gas is in a hydrostatic balance. 
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p dr dr 

where P and p arc the gas pressure and density, respec- 
tively, and the gravitational potential and A'/dm('') is 
the dark matter mass contained in a radius r. We solve 
equation (1) with an equation of state as P oc p'^ to ob- 
tain the initial gas distribution, where 7 = 1-1-1 /N for 
the polytrope index N. We consider the NFW distribu- 
tion and approximated King distribution (referred to as 
King distribution hereafter) for the dark matter potential 
as shown in figure 1. 

The gas temperature profile thus obtained for the NFW 
potential can be expressed in a form 
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where rg is the characteristic radius of the NFW distribu- 
tions expressed as 
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and dcPc is a constant related with the dark matter mass 
Mdm through 
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Fig. 1. Dark matter density profiles with = 400 Icpc for 
the NFW (solid) and Vc = 100 kpc for the King (broken) dis- 
tribution taken into calculations here (see text). The King 
profile with rc = 400 kpc is also shown for reference. The 
three curves are normalized so as to pass the same point at 
r = 400 kpc. 

The gas density profile is given by 
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For the King potential, similarly we obtain the gas tem- 
perature as 
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Here Tc is the core radius of the King distribution ex- 
pressed as 
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where po is the dark matter density at the center. The 
dark matter mass is given by 
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The gas density profile for the King potential is given by 
equation (6) as well. 

The density profiles, p{r)/p{0), calculated from the 
above equations are shown by the broken lines (denoted 
by NCC) in figure 2 with r^ = 400 kpc for the NFW and 
Tc = 100 kpc for the King dark matter potentials. These 
radii are chosen for the /3-model core radius of surface 
brightness (cx p^T^/^) to give nearly the same radii of 
100-200 kpc as each other. The equations (2)-(6) are 
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equivalent to those (12) and (42)-(46) in Suto, Sasaki & 
Makino (1998), and therefore the broken lines in the left 
panel (NFW NCC) of figure 2 arc the same results as 
theirs. We calculate the temperature profile, T{r)/T{0), 
as well for calculations in the following section. 

2.2. Quasi-hydrostatic gas distribution under radiative 
cooling 

Now we consider gas distribution of a cool core (CC) 
that significantly undergoes radiative cooling. The gas is 
then no longer in the hydrostatic balance but should be 
described by hydrodynamics, since inflow toward the inner 
region where the cooling rate is higher would occur so as 
to compensate the pressure decrease due to radiation loss 
(Fabian 1994). Hence, the gas temperature and density 
profiles of the cluster core deviate from those polytropic 
profiles given in the last subsection. 

In the early cooling phase, however, inflow is yet mild 
and quasi-hydrostatic condition is attained marginally un- 
til the initial cooling time TcooI- Here TcooI is defined with 
the density and temperature at the center of a virialized 
cluster before cooling and is typically ~ Gyr (Akahori & 
Masai 2006). We consider that the pressure decrease at 
a radius due to radiation is compensated immediately by 
gas inflow from the adjacent outer (and hotter) region, or 
the local inflow rate is controlled by the local cooling rate, 
and thus quasi-hydrostatic condition is attained. 

The profile of the cooling gas temperature, T'{r), is then 
expressed approximately in a form (Masai & Kitayama 
2004) 



d\uT\ 9r 1 lo 
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at < r < 

^cooij when variation of T(r^ (NCC) is enough 
small in this range; it should be noted that this equa- 
tion was obtained originally for isothermal mass distribu- 
tion of the virial temperature Tq. Here, TcooI is a cooling 
core radius in which the local cooling time is shorter than 
the Hubble time, as tcooi(r) = 3A:T(r)/[2p(r)A] < H^^ at 
T < TcooI and icooKj^cooi) = ^ where pA represents the 
radiative cooling rate per unit time. We assume [5t ^l, 
i.e., the gas and dark matter have the same tempera- 
tures as each other. In observations, the values of /? 
implied by the density profile fall around 0.6 less than 
unity. This may, however, be a result inherent to the 
/3-model because it is normally applied to observed pro- 
files within a radius significantly smaller than the viral 
radius (see Akahori & Masai 2005). We put T(0) = To 
for normalization and calculate T'(r)/T(0) with a bound- 
ary condition T'(rcooi) = T(rcooi)- It should be noted that 
T{rcoo\)/T{0) depends on the potential profile and 7, and 
thus T'(r)/T(0) does as well, though not explicitly seen 
in equation (11). 

Once the gas inflow takes place to be steady and the 
quasi-hydrostatic condition is attained, we may have 
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lively, of the cool core. The ram pressure is negligible 
compared to the thermal one in this balance. Thus, for 
given T'(r) the density profile 



d\np' 
dlnr 



dlnT' 
dlnr 



T 
T 



dlnp 
''^dlnr 



(13) 



follows; note that 7 is defined for P (or T) and p. In the 
cool core concerned, T/T' varies with r from 1 at TcooI 
to ~ 3 at r ^ 0, and the gas is kept (quasi-) hydrostatic 
through the core. Actually, for infiow to occur, the left- 
hand side must be once larger than the right-hand side in 
equation (12); note that dP/dr is negative. 

Using equation (13) with T(r) and p(r) obtained in the 
last section and T'{r) from equation (11), we calculate 
p'{r)/p{Q) with p'{rcoo\) = pircooi), and show in figure 2 
by the solid lines (denoted by CC) for a typical cooling 
radius of TcooI = 100 kpc (e.g., Peterson & Fabian 2006). 
The initial gas density is taken to be a polytropic distribu- 
tion of 7= 1.001 or 1.66 (see Section 2.1) and represented 
by the broken line (NCC). The virialized gas of 7~ 1 poly- 
trope may be attained, e.g., if merger shocks are radiative 
or thermal conduction works through the core in the for- 
mation history of the cluster. Note that quasi-hydrostatic 
cooling (Masai & Kitayama 2004) predicts a steeper pro- 
file than the hydrostatic one of 7 = 1.001 polytrope. 

In figure 2, the dotted lines (denoted by (3) represent 
the best-fit /3-model, n(r) = no [1 -I- (r/rcgas)^]"^*^^^ with 
n and r^gas being the density and the core radius of gas, 
to mimic the density of CC at 5 < r/kpc < 100. The 
parameters are given in table 1, where also is given the 
parameters for the initial non-cool cores in the range 5 < 
r/kpc < 400. At r < rcooi the /3-model can reproduce the 
density of CC fairly well for the NFW and quite well for 
the King potential; systematically, the /3-model gives a flat 
profile toward the center than calculated for the NFW. 
This result may suggest that analysis with the /3-model 
is useful also for cool cores if the gas is quasi-hydrostatic, 
though not relevant to its original meaning. 

Figure 3 shows the pressure profile, p{r)T{r)/[p{0)T{0)] 
(broken lines, NCC) and p'(r)T'(r)/[/9(0)T(0)] (solid lines, 
CC). It is seen that in the quasi- hydrostatic cooling core 



(r 



< 



100 kpc) the pressure profile is steeper for 



where p' and P' are the gas density and pressure, respec- 



the initially isothermal (7 = 1.001) gas than the 7 — 
1.66 gas, and for the NFW potential than the King 
one well inside TcooI- Also shown by the dotted line is 
the best-fit generalized NFW (GNFW, hereafter; Nagai 
et al. 2007) pressure obtained for the REXCESS sam- 
ple of clusters (Arnaud et al. 2010 Eqs. 11-12), p{x) = 
/'o[(c50ox)-«(l + (c5ooa;))^)]-("-'')/^ for (Po,C5oo,0,C,^) = 
(8.403/lTO/^ 1.177,0.3081,1.0510,5.4905). The GNFW 
profile is normalized so as to pass the intersection of 
7 = 1.001 and 7 = 1.66 curves of quasi-hydrostatic cooling 
gas (CC). 

Figures 4 shows the entropy profile, 
T(r)/p2/3(r)/[r(0)/p(0)2/3] (broken lines, NCC) and 
r'(r)/p'2/3(r)/[T(0)//9(0)2/3] (solid lines, CC). While the 
initial polytropic gas (NCC) exhibits a fiat or shallower 
profile, the cooling gas (CC) exhibits entropy roughly 
cx r at 0.1 < r/rcooi < 1 for both the NFW and King 
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Fig. 2. Gas density profile of quasi-fiydrostatic cool cores, p'{r)/p(0) (solid lines, CC) with r^^^l = 100 kpc and the initial non-cool 
ones, p{r)/p{0) (broken lines, NCC) of 7 = 1.66 (thick lines) or 1.001 (thin lines) polytrope. The dotted lines represent the best-fit 
/3-model (see text). Particularly in the right panel, the best-fit /3-modcls (thin and thick dotted lines) reproduce the calculated CC 
density profiles (thin and thick solid lines) well and almost perfectly overlap with them for both 7 = 1.001 and 1.66, respectively. 




Fig. 3. Gas pressure profile of quasi-hydrostatic cool cores, p' {r)T' {r)/[p{0)T{0)] (solid lines, CC) with rcool = 100 kpc and the 
initial non-cool ones, p(r)T(r)/[p(0)T(0)](broken lines, NCC) of 7 = 1.66 (thick lines) or 1.001 (thin lines) polytrope. The dotted 
lines represent the GNFW pressure profile obtained for REXCESS sample by Arnaud et al. (2010) (see text). 
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Fig. 4. Entropy profile of quasi-hydrostatic cool cores, T'(r)/p'2/3(r)/[T(0)/p2/^(0)] (solid lines, CC) with r^aol = 100 kpc and the 
initial non-cool ones, T(r)/p2/3 (r)/[T(0)/p2/^ (O)](broken lines, NCC) of 7 = 1.66 (thick lines) or 1.001 (thin lines) polytrope. The 
upper and lower dotted curves in each panel indicate the best-fit entropy profiles for the ACCEPT sample with low {Kq < 50 keVcm^) 
and high {Kq > 50 keVcm^) central entropy values (Cavagnolo ct al. 2009), respectively. 



Table 1. /3-model parameters for density profiles of non-cool and cool cores. 







7 


= 1.001 




7 


= 1.66 




Potential 




no 


'"cgas [kpc] 




no 


^cgas [kpc] 


13 


NFW 


NCC 


0.89 (1.00) 


64 


0.33 


0.95 (1.00) 


104 


0.25 


NFW 


CC 


5.41 (6.24) 


25 


0.55 


4.05 (4.43) 


33 


0.45 


King 


NCC 


1.00 (1.00) 


100 


0.53 


0.99 (1.00) 


164 


0.44 


King 


CC 


5.02 (5.05) 


51 


1.0 


4.10 (4.13) 


58 


0.81 



Central density normalized to that for the non-cool core. The value in parenthesis denotes the central surface 
brightness obtained from the present calculation. 



potentials. The best-fit profiles obtained by Cavagnolo et 
al. (2009) for the ACCEPT cluster sample are also shown 
with the dotted lines so as for their high central entropy 
Kq (> 50 kcV cm^) to be unity at r = 0. Note that their 
entropy profile consists of a power-law plus a constant, 
K{t) = Kq + iirioo(r/100 kpc)", and the best-fit parame- 
ters are {Ko,Kioo,a) = (16.1, 150, 1.20), (156, 107, 1.23) for 
clusters with Kq < 50 and Kq > 50 keVcm^, respectively. 

Comparison of the calculated density, pressure and en- 
tropy with observations is discussed in the next section. 
It should be noted here that TcooI is a point of discontinu- 
ity since we concern cooling within rcooi leaving the initial 
polytropic value as is outside. 

The picture discussed in this section is valid as long 
as sound crossing time is enough shorter than TcooI- 
In other words, without any heating process, sound 
speed decreases with cooling and eventually becomes un- 
able to readjust pressure against the gravitational force. 
Flattening of p'{r) toward the center is due to (quasi-) 
hydrostatic balance, and actual p'{r) could be steeper in 
hydrodynamics. Note that T' or p' here is not applica- 
ble to the deep core where the gas accumulates, and also 



recall that we ignore the gas mass for the gravitational 
potential. In the quasi-hydrostatic cooling phase, equa- 
tion (12) or (13) can be a good approximation for typical 
cool cores. The results here are useful for /3-model based 
or such hydrostatic balance based analysis, as discussed 
in the following section. 

3. Discussion 

3.1. Comparison with observed core radii and density 
profiles 

As seen in figure 2 or table 1, the core radius in the 
quasi-hydrostatic cooling phase is about one-third of the 
initial radius of polytropic gas in the virial equilibrium. If 
clusters which are relaxed or virialized but not yet cooling- 
affected have ^ 200 kpc cores, the clusters in the quasi- 
hydrostatic cooling stage appear to be those of ^ 60 kpc 
cores. This result is consistent with the hydrodynamical 
calculations by Akahori & Masai (2006), who analyzed the 
calculated clusters with the /3-model. 

In figure 5, the gas density profiles are compared with 
the best-fit /3-model profiles for clusters analyzed in Ota 
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Fig. 5. Comparison of calculated gas density profiles with observations. The solid and broken black lines represent the density 
profiles of quasi- hydrostatic cool cores (CC) and the initial non-cool cores (NCC), respectively, for the NFW (left) and King (right) 
potential. The calculated density profiles are normalized such that the cooling time at rcool is equal to 10 Gyr. The blue and red 
lines show the electron density profiles of clusters with small (rc < 100 kpc) and large (>100 kpc) core radii, respectively (taken 
from the /3-model analysis in Ota et al. 2006.) 



et al. (2006). Note that the systematic error of the cen- 
tral gas density derived under the isothermal /3-model is 
estimated to be small (Appendix 1). which does not affect 
the present discussion. In figure 5, the solid (CC) and 
broken (NCC) black lines represent the density p' of the 
quasi-hydrostatic cool core and p of the initial polytropic 
non-cool core, respectively. These calculated density pro- 
files are normalized such that the cooling time at TcooI is 
equal to 10 Gyr in the figure^. 

One can see that the initial density p (broken lines) falls 
around the upper bounds of the densities of the larger-core 
(rc > 100 kpc) clusters, likely most relaxed but not yet 
cool-core ones. The density p' (sohd lines) of the quasi- 
hydrostatic cool core, on the other hand, lies between 
the mid- and high-values of the densities of small-core 
(rc < 100 kpc) clusters. We also found from literatures on 
spatially-resolved X-ray spectroscopy of clusters (Allen et 
al. 2001; Bonamente et al. 2006; Zhang et al. 2007; Zhang 
et al. 2008; Snowden et al. 2008) that 11 out of 13 small- 
core clusters plotted in figure 5 have a clear temperature 
decrease at the center while none of 17 large-core clusters 
has such a temperature drop. 

The above facts may support that the density increase 
in the cool cores can be accounted for by quasi-hydrostatic 
cooling and our calculation gives a typical density of cool 
cores. However, there exists a large scatter and some 
observed clusters have central densities higher than pre- 
dicted. In the context concerned here, the scatter is as- 

^ The temperature at the cooling radius is assumed to be 
fcT(r(.ooi ) = kT' (r^ool ) = 5 keV in the calculation of cooling time, 
whose dependence on fcT is, however, sufficiently small that it 
does not affect the present discussion. 



cribed to that in the potential profile and/or the initial 
infiow condition and time elapsed. 

3.2. Comparison with observed pressure and entropy pro- 
files 

Arnaud et al. (2010) discussed universal pressure pro- 
file for the REXCESS cluster sample, and obtained the 
best-fit profile (their equations 11-12) based on GNFW by 
Nagai et al. (2007). It is seen from figure 3 that the best- 
fit GNFW profile is close to the overall (10 < r/kpc < 200) 
slope of quasi-hydrostatic cooling gas of initially 7 = 1.001 
polytrope under the NFW potential. Compared to the 
GNFW, a depression in the mid-part could arise from dis- 
continuity in p' and T' across TcooI (see §2.2). On the other 
hand, the difference at ^ 200 kpc is due likely to that in 
the dark matter potential or rg, or cluster gas is not sim- 
ply hydrostatic otherwise. In any case of 7, it is unlikely 
that hydrostatic gas like initial polytrope here (NCC) can 
account for such a steep profile as the GNFW in the inner 
core. 

Cavagnolo et al. (2009) derived radial entropy pro- 
files of ICM for 239 clusters with the Chandra data (the 
ACCEPT sample) to find that most entropy profiles are 
well fitted by the model consisting of a power-law plus a 
constant K{r) = Kq + Kioo{r /lOO kpc)" (see also §2.2). 
They also showed that the distribution of central entropy 
Kq is bimodal, which peaks at Kq ~ 15 keVcm^ and ~ 
150 keVcm^. From a comparison of the calculated entropy 
lies with the best-fit models for the ACCEPT sample (fig- 
ure 4), we find that the model for clusters with higher 
central entropy (A'o > 50 keVcm^) is close to the initial 
entropy profiles for 7 ~ 1.001. As for A'o < 50 keVcm^, 
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there is a reasonably good agreement with the CC cases 
in the central decrement, which is by about an order of 
magnitude lower than that of the initial NCC profiles. 

Moreover, for the scaled temperature and the density, 
Arnaud et al. (2010) found that their deviations from the 
average scaled profile are anti-correlated with each other 
in the core, r/i?5oo < 0.2 (figure 3 in their paper); the anti- 
correlation is more clearly seen for cool core clusters. This 
behavior is in favor of quasi-hydrostatic cooling picture. 

3.3. /3 -model fit to surface brightness profiles: physical 
basis of applying double-j3 model to cool core clusters 

The density profile of the quasi-hydrostatic cool core 
can be mimicked by the /3-model with a core radius smaller 
than that for the initial non-cool core, as demonstrated in 
figure 2. This suggests that a cool-core cluster can be still 
fitted to either of the /3-model or double /3-model, depend- 
ing also on the available range of radius. Observationally, 
the double-/? model is frequently used to fit the X-ray data 
of clusters having compact cores or cool cores (e.g., Ota & 
Mitsuda 2004; Santos et al. 2010). In this subsection, we 
examine whether X-ray surface brightness expected under 
quasi-hydrostatic cooling can be fitted by the conventional 
single-component or two-component /3-models. 

We calculate the surface brightness of free- free emission 
(oc p'^T'^/^ or oc p^T^I"^ of cool or non-cool core). As 
readily seen in figure 6, that the gas core appears smaller 
in the quasi-hydrostatic cool core than the initial non-cool 
one. We performed the /3-model fitting to the calculated 
surface brightness profiles with the form S[r) = S'o[l + 
{rlrfi^YY'^^^^l'^ . The radial range of 5 < r/kpc < 400 is 
used. The best-fit parameters are given in table 2. The 
surface brightness is enhanced by a factor of 2.4-6.5 and 
the core radius becomes smaller by a factor of ^ 2.5-10 
for the quasi-hydrostatic cool core in comparison with the 
initial non-cool core. In the lower panels of figure 6, the 
ratio of the calculated surface brightness profiles to the 
best-fit /? models are shown. The single-component /?- 
model gives a reasonably good fit to the NCC case, while 
it leaves a systematic deviation from that of the CC case. 
The chi-squared value defined as = ~ '5'(^)]^ is 

also given in table 2. Here l(r) represents the calculated 
surface brightness profile and ^ sums up all data points. 

Next, the double-/3-model, 5(r) = Y!i=\ + 
{rlrf,t,if\~''^^'^^l'^ , is applied to the CC profiles. In this 
fitting, the inner and outer slope parameters are tied (i.e., 
/?! =1^2) because of a large uncertainty of the inner slope 
and allowed to vary within 0-2. The two core radii, rst,! 
and rfit,2 are restricted to a range of 5-400 kpc. As is 
seen from figure 7, the surface brightness profile of clus- 
ters under radiative cooling can be well approximated by 
this two-component model. In fact, the systematic devia- 
tion in the single-^ model fitting is significantly reduced 
and the values become smaller in this case (table 2). 

The above result gives a physical basis for the use of 
double-/? model in the analysis of X-ray surface bright- 
ness of cool core clusters. We should note that the result- 
ing rfit^i or rfit.2 in this analysis are not necessary equal 



to the core radius assumed for the potential distribution, 
Tc, or that expected through the relationship, 7'c ~ Q.22rs 
(Makino et al. 1998). 

The inner component with small core size of r^t,! ^ 
20 — 70 kpc (table 2) can be ascribed to radiatively cool- 
ing gas at the center. The present result will also explain 
the fiducial emission measure profile obtained by Arnaud 
et al. (2002) as well as the double-peaked core-size dis- 
tribution discovered by Ota & Mitsuda (2002); Ota & 
Mitsuda (2004), whose peak values are r^ = 50 kpc and 
200 kpc. 2 

3.4^. Implication on the time scale of ICM thermal evo- 
lution 

The fraction of cool core clusters will give another im- 
portant clue to understand the formation of cool cores 
since it should refiect any timescale relevant to their evo- 
lution. 

Akahori & Masai (2006) calculated the probability dis- 
tribution by estimating the time during which a cluster 
would have a certain core size, to confirm the double- 
peaked nature of core radius. They suggest that the gas 
rapidly cools after the quasi-hydrostatic condition breaks. 
In other words, relaxed clusters with small core radii 
currently observed are more or less those in the quasi- 
hydrostatic cooling stage, which terminates at ~ TcooI- 
Since TcooI is estimated to be smaller than the Hubble 
time H^^ ^ 10 Gyr for a typical cluster, some heating 
is needed to sustain the system otherwise it would be 
unseen in ^ Gyr after virialized. Practically, however, 
heating due to mergers is likely invoked in the cluster 
evolution, as suggested also from a point of view of ra- 
dio halos (e.g., Rossetti et al. 2010). The clusters of core 
radii > 400 kpc in the core-size distribution arc attributed 
possibly to mergers from their morphology; it should be 
noted that quasi-hydrostatic cooling is a picture for a re- 
laxed gravitational system. Also other heating/pressure 
sources, e.g. by black hole activity, may work. In fact, 
many of so-called cooling flow clusters have central cD 
galaxies hosting black holes. Radio mini-halos observed 
from these clusters may suggest the considerable contri- 
bution of the black holes. 

Utilizing the technique of the X-ray fundamental plane, 
Ota et al. (2006) investigated the evolution of X-ray ob- 
servables for the distant cluster sample. They showed that 
the principle axis of the plane is parallel to the tcooi-axis, 
suggesting the cooling time is a parameter to control the 
core structure. They also noted that the rate of radiative 
energy loss is likely to be kept nearly constant for a signifi- 
cant duration of time even after the onset of cooling in the 
relaxed systems. This indicates that some steady state is 
attained for the gas of many small-core clusters. Those 
results are consistently understood within the framework 

^ Later, a similar distribution has been shown independently by 
Hudson et al. (2010) using the high-resolution Chandra data on 
a nearby flux-limited sample. Notice that the detailed statistics 
of the cool-core formation seems to admit of further investiga- 
tion (e.g., Eckert et al. 2011 and references therein), which is, 
however, beyond the scope of the present paper. 
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Fig. 6. Surface brightness of free-free emission calculated for quasi-hydrostatic cool cores (solid lines, CC) and the initial non-cool 
cores (broken lines, NCC). In the upper panels, the best- fit single /3-model are also indicated with the red/blue dotted curves for 
the NCC/CC cases with 7 = 1.001 and the red/blue dash-dotted curves for the NCC/CC cases with 7 = 1.66, respectively. In the 
lower panels, the ratio of the calculated surface brightness profiles to the best-fit models are plotted. 



Table 2. /3-model parameters for surface brightness profiles in the non-cool and cool cores. 



Potential 
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rfit,i [kpc] 
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''fit, 2 [kpc] 
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single-/? 
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double-/? 


5.14 
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0.37 








0.41 
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0.49 


King 


1.66 
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single-/? 


3.42 (3.48) 
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17.45 


King 


1.66 
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2.74 


74 


1.19 


0.69 


336 


1.19 


0.30 



Central surface brightness normalized to that for the non-cool core. The value in parenthesis denotes the central 
surface brightness obtained from the present calculation. 
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Fig. 7. Surface brightness of free-free emission calculated for quasi-hydrostatic cool cores (solid lines, CC). In the upper panels, 
the best-fit single-/? and double-/? models are indicated with the blue dotted and magenta dashed curves. The contributions of two 
components of the double-/? model are also shown with the red and blue curves, respectively. In the lower panels, the ratios of 
the calculated surface brightness for the cool core to the best-fit single-/? and double-/? model are shown with the blue dotted and 
magenta dashed curves, respectively. 
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of quasi-hydrostatic model. 

Finally, it should be noted that the present calculation 
shows the evolution of density profile in the cluster core 
when the initial NCC profile is given by p{r)/p{0) in §2.1: 
a different initial condition may result in a different CC 
profile. We ignore the mass of inflow gas through our 
calculation. In the very late stage of quasi-hydrostatic 
cooling, however, the inflow gas could alter the gravita- 
tional potential in the inner core region to a steeper pro- 
file from the initial equilibrium one. This would break the 
quasi-hydrostatic condition and lead to the global cool- 
ing flow eventually, unless some counter pressure works 
against, e.g., by heating of the gas. Further investigation 
is needed to obtain the total view on the thermal evolu- 
tion of cluster core, including the possibility of heat input 
due to mergers as mentioned above. 

4. Summary 

1. We have calculated the density profile of cool core 
of intracluster gas in the quasi-hydrostatic cooling 
phase (Masai & Kitayama 2004) assuming that the 
gas is initially in the virial equilibrium state within 
the dark matter potential represented by the NFW 
or King distribution and has a polytropic profile 
with 7== 1.001 or 1.66. 

2. In the quasi- hydrostatic cooling phase, the gas den- 
sity increases by a factor of 4-6 at the cluster center 
while the temperature decreases to be about one- 
third in comparison with their initial polytropic pro- 
files (non-cool core, NCC). Accordingly, the core ra- 
dius of the cool core (CC) cluster appears more com- 
pact than the initial polytropic one. This result is 
consistent with the hydrodynamical simulations by 
Akahori & Masai (2006). 

3. From the comparison with X-ray observations, it is 
found that the density profile of CC falls between the 
mid- and high-values of small-core {r^ < 100 kpc) 
when normalized so that the initial density profile 
has the cooling time equal to 10 Gyr at the cooling 
radius. 

4. The pressure and entropy profiles are derived for 
the quasi- hydrostatic cooling gas and their CC/NCC 
profiles are compared with the generalized NFW 
pressure profile for the REXCESS cluster sample 
(Arnaud et al. 2010) and the best-fit entropy pro- 
file models for the ACCEPT sample (Cavagnolo et 
al. 2009), respectively. 

5. The X-ray brightness profiles are calculated for clus- 
ters having either the NFW or King potential as- 
suming the free-free emission. The resulting CC 
profiles are found to be well represented by the con- 
ventional double /3-model including the central emis- 
sion with a compact core size. This gives a physical 
basis of applying the double-/? model to observed 
CC clusters. The inner component of the double- 
/3 model, which can be ascribed to the radiatively 
cooling gas at the center, will also explain the fidu- 
cial emission measure profiles (Arnaud et al. 2002) 



as well as the double-peaked core-size distribution 
(Ota k Mitsuda 2002; Ota & Mitsuda 2004). 
6. Implication on the time scale of ICM thermal evolu- 
tion is discussed. Since the quasi-hydrostatic stage 
will terminates at TcooI ^ Gyr < Hq^ for a typi- 
cal cluster, some heating is needed to sustain the 
system. The mass of inflow gas is ignored in the 
present calculation, however, it will alter the the 
potential distribution of the inner core to a even 
steeper one in the late stage of cooling, which would 
break the quasi-hydrostatic balance and lead to the 
global cooling flow. Further investigation is needed 
to construct the total view on the thermal evolu- 
tion, including the possibility of heat input due to 
mergers. 

This work was supported in part by the Grant-in- Aid 
for Young Scientists by MEXT, No. 22740124 (NO) and 
by the Grant-in- Aid for Scientific Research 18540241 from 
JSPS (KM). 

Appendix 1. Systematic error of central gas den- 
sity 

Ota et al. (2006) derived the electron density by the 
isothermal /3-model fitting to the surface brightness pro- 
files. Since the temperature dependence of X-ray surface 
brightness is generally small (cxT^/^), the systematic error 
of central gas density coming from the temperature drop 
at the cluster center is expected to be small. Actually, 
wc made a comparison with the central gas density ob- 
tained from the XMM-Newton observations having im- 
proved resolutions (Zhang et al. 2007) and confirmed that 
there is a good agreement between the two samples as 
shown in figure 8. The figure includes 7 cool core clusters 
as well as 2 non-cool core clusters. Thus the evolution of 
gas density profile can be further inferred by comparing 
those calculated for the cool cores and non-cool ones in 
the present work with the /3-modcl profiles presented in 
Ota et al. (2006). 
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